quadrature subroutine

平方子程序